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Abstract 

We study the integer sequence v n of numbers of lines in hypersurfaces of degree 
2n — 3 of P n , n > 1. We prove a number of congruence properties of these numbers 
of several different types. Furthermore, the asymptotics of the v n are described 
(in an appendix by Don Zagier). An attempt is made at a similar analysis of 
two other enumerative sequences: the numbers of rational plane curves and the 
^ . numbers of instantons in the quintic threefold. 

We study the sequence of numbers of lines in a hypersurface of degree D = 2n — 3 of 
^ P™, n > 1. The sequence is defined by (see e.g. [B]) 

v„:= c 2ri - 2 (Sym D (5), (1) 

1 JG(2,n+l) 

where G(2,n + 1) is the Grassmannian of C 2 subspaces of C n+1 (i.e. projective lines in 
P n ) of dimension 2(n + 1 — 2) = 2n — 2, Q is the bundle of linear forms on the line (of 
>- . rank r = 2, corresponding to a particular point of the Grassmannian), and Sym° is its 

Dth symmetric product - of rank (^^I^ 1 ) = D — 1 = 2n — 2. The top Chern class (Euler 
class) C2n_2 is the class dual to the 0-chain (i.e. points) corresponding to the zeros of 
the bundle Sym°(Q), i.e. to the vanishing of a degree D equation in P n ; this is the 
geometric requirement that the lines lie in a hypersurface. 
The integral (P) can actually be written as a sum: 



E 



Ua=o( aw i + (D- a)wj) 



(2) 



where wq, . . . ,w n are arbitrary complex variables. This is a consequence of a localisation 
formula due to Atiyah and Bott from equivariant cohomology, which says that only the 
(isolated) fixed points of the (C*) n+1 action contribute to the defining integral of v n . 
Hence the sum. For the first few values of n computation yields t>2 = 1, v% = 27, 
v 4 = 2875, v 5 = 698005, v 6 = 305093061. 

D. Zagier gave a simple proof that the right hand side of (J2J) is independent of 
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Wo, . . . ,w n (as it must be for (J2j) to hold), and that in fact it can be replaced by the 
much simpler formula 

2n-3 

[l-x)T[(2n-3-j + jxj\ (3) 

J X n ~ 

3=0 

where the notation [. . ,] x n means the coefficient of x n . In fact, formula (J2J) was proved 
in a very different way using methods from Schubert calculus by B.L. van der Waerden, 
who established it in part 2 of his celebrated 20 part 'Zur algebraischen Geometric' 
series of papers |T9"1 |2*U] . The number of linear subspaces of dimension k contained in 
a generic hypersurface of degree d in P n , when it is finite, can be likewise expressed as 
the coefficient of a monomial in a certain polynomial in several variables, see e.g. [T5J 
Theorem 3.5.18]. 

Zagier also gave the formula 

[27, x2 „ 7/ ,/ 9 111 9999 \ 

where the right hand side is an asymptotic expansion in powers of n _1 with rational 
coefficients that can be explicitly computed. The proof of this formula, as well as the 
derivation of (JHJ) from (j2J, can be found in the appendix. 

The remaining results, summarized in Theorem^ and Theorem are concerned with 
congruence properties of the numbers v n . In this context it turns out to be convenient 
to define v\ = 1 (even though there is no such thing as a hypersurface in P 1 of degree 
— 1) and even more remarkably Vq = —1. We do not doubt that the congruence results 
presented here form only the tip of an iceberg. 

A first version of this paper was single authored by the first author. Indeed, the 
present version of this paper is similar to the first one, except for sections 2 and 3 which 
have been greatly revised and expanded by the second author. The conjectures outside 
these two sections are due to the first author alone. Sections 4 and 5 were revised by 
both the second author and Don Zagier. 



1 Introduction 

The motivating idea behind this paper is the expectation that certain problems in enu- 
merative geometry are coupled to modularity. This is a recurrent theme in string theory, 
where partition functions have often an enumerative interpretation as counting objects 
(instantons, etc) and must satisfy the condition of modularity covariance in order to 
obtain the same amplitude when two worldsheets have the same intrinsic geometry. 

Modular forms, as is well known, have Fourier coefficients satisfying many interest- 
ing congruences (think of Ramanujan's congruences for partitions, or for his function 
r(n)). The same can happen for the coefficients of expansions related to modular forms, 
e.g. the expansions y = A n x n obtained by writing a modular form y (locally) as a 
power series in a modular function x. For instance, the famous Apery numbers related 
to Apery's proof of the irrationality of £(3) are obtainable in this way [2j and satisfy 
many interesting congruences |18j . The numbers appearing in the context of mirror 
symmetry, Picard-Fuchs equations for Calabi-Yau manifolds, Gromov-Witten invariants 
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and similar problems of enumerative geometry are sometimes related to modular forms 
and sometimes not, so we can reasonably hope for interesting congruence properties in 
these contexts also. In Section 2 we shall find astonishingly many congruences for our 
sequence v n . We shall first draw a few tables for congruences mod 2,3,4,5 or 11, and 
then summarize the observed congruences. In Section 3 we prove those congruences by 
elementary means starting from (J3j), and a few conjectures will be formulated. Sequences 
of numbers coming from modular forms also often have interesting asymptotic properties 
and we therefore wish to study this, too. In Section 4 we find the asymptotic properties 
of the v n numerically by using a clever empirical trick shown to us by Don Zagier which 
we call the asymp^ trick. (A rigorous proof of these asymptotics, as already mentioned 
above, was also provided by him and is reproduced in the appendix.) Section 5 presents 
congruences and asymptotics for two further examples of enumerative sequences, with- 
out proofs. The sequence of rational curves on the plane and the sequence of instantons 
on the quintic threefolds partly partly mimic the behaviour of the original sequence of 
lines in hypersurfaces. 



2 Congruences 



We will consider the sequences {v n (mod k)}^ =1 for some small values of k\ that is, we 
study the reduction of the integers v n modulo k. It turns out to be instructive to order 



the v n modulo k in a table. Each table has k columns. The ith. column {i 
gives the values of vik+i (mod k) {I > 0). 

For instance, the first few tables at k = 2, 3, 4, . . . look like 

11222000 

k = 2: ; ; ; ; ; ; ; ■ ■ ■ k = 3-. 1 1222000 
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k) 
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The first table (k = 2) says that all the v n are odd integers. We shall mostly be inter- 
ested in the tables for prime k. Here is a typical prime table: 



k = 11: 
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Study of these and other tables led us to formulate a number of conjectures, most of 
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which we were able to prove. An overview of these results is given in Theorem 21 

Theorem 1 The following holds for the tables of v n mod k: 

1. All v n are odd. 

2. The first two rows of each table are equal. 

3. If k is even, then rows k/2 + 1 and k/2 + 2 are equal. 

4. For k odd, row (k + 3)/2 contains only zeros. 

5. For k prime, the first two rows start with 1,1 followed by k occurrences of k — 1. 

6. For k > 2 prime, the last (k — l)/2 slots of the first column vanish. 

7. For k > 2 prime, there is a block of zeros at the bottom (after (k — l)/2 columns), of 
height (k — l)/2 and width (k + 3)/2. 

8. For k = 2 q , all rows are constant and in two fold way sweep out all odd residues, i.e. 
for every odd integer a with 1 < a < 2 q there are precisly two rows that have only a as 
entry. 

9. For k = 2 q > 2 the entries in the rows 1, 2, 2 9_1 , 2 q ~ l + 1,2 9_1 + 2, 2 q - 1 equal, 
respectively, 1, 1, 2 q ~ l - 1, 2"- 1 + 1, 2 q ~ l + 1, 2« - 1. 

10. For k = 2 q > 2 the entries in row a and row a + 2 9 " 1 differ by 2 q ~ 1 (mod 2 9 ). 

Proof. These ten claims are proved in, respectively, lemmas El El El HI HI & El HH Ed EDI 
&ESlI23and[221 □ 



On computing the reductions of t>i, . . . ,t> 32 modulo 32 one finds by part 8 of this 
theorem that for k = 4 the table has constant rows 1,1,3,3, for k = 8 constant 
rows 1,1,3,3,5,5,7,7, for k = 16 constant rows 1,1,11,11,5,5,7,7,9,9,3,3,13,13, 
15, 15 and for k = 32 constant rows 1, 1, 27, 27, 21, 5, 7, 23, 9, 9, 19, 19, 29, 13, 31, 15, 17, 17, 
11, 11, 5, 21, 23, 7, 25, 25, 3, 3, 13, 29, 15, 31. Thus, for modulus 2 q with q < 3 we observe 
that pairs of values occur and that these, moreover, are in ascending order. For q = 4 
the values still come in pairs, but the order is no longer ascending. For n > 5 it turns 
out that pairs with equal values become sparser and sparser. Notice that in the above 
cases for every modulus all odd values are assumed exactly twice. By part 8 this always 
happens. Thus, given an odd integer a and any integer q > 1 there are infinitely many 
integers m such that v m = a(mod 2 q ) (or put more succinctly: modulo powers of two 
the sequence v n is equidistributed over the odd residue classes). 

For k prime, often v n = 0(mod k) for trivial reasons. It then makes sense to consider 
divisibility of v n by higher powers of k. Our deepest result in this direction is provided 
by the following theorem. 

Theorem 2 1. If p > 5 is a prime, then 

V p+3 = -2p 3 (mod p 4 ) and Vp+a = 2p 3 (l -p)(p- 1)!4 P_1 (mod p 5 ). 

2. Let r > 1 and p > 2r + 1 be a prime. Then 

v P _+i +rp = C rP 2r+2 (modp 2r+3 ), (5) 

where 



j=0 

the integers bj r are defined implicitly by Y\ 2 a=i(^ r + 1 — a + ax) = Y^[=Q°j , and 
((«)). :=IE=i(« + 2j -2). 
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Remark. Note that &2r-j,r = bj :T . Numerical experimentation suggests that the numera- 
tor of c r always equals a power of 2 and that the congruence (jSJ holds for all odd primes. 
Below we record some values of c r . 



1 2 3 4 5 

c r- _oi 103125 210171535 1308348857025 11660783598520749 

° 8 64 1024 16384 

3 Proofs of the theorems 
3.1 Some generalities 

First recall from the elementary theory of finite fields of order p, that 

p-i 

a; p_1 — 1 = Y\( x ~~ 3) ( m °d p)- 

3=1 

(Here and below the letter x denotes a variable.) By substituting x = one obtains 
Wilson's theorem: 

(p — 1)! = —1 (mod p). 

We also recall the freshman's identity (a + b) p = aP + IP (mod p), from which we infer 
that if f(x) e Z[x], then f(x) p = f(x p ) (mod p). These results will be freely used in the 
sequel, without further referring to them. 

Lemma 1 We have v n = (mod (2n — 3) 2 ). 



Proof. The term with j = in Q equals 2n — 3. The term with j = 2n — 3 equals 
(2n — 3)x. Hence v n = (2n — 3) 2 w n , where w n = (1 — x) Y\^=\ i {^ n ~ 3 — j + jx 



. □ 

n-2 



The following result was first noticed by D. Kerner. An alternative, slightly longer, 
proof was given by M. Vlasenko. 

Lemma 2 We have v n = (mod (2n — 3) 3 ). 

Proof. It suffices to show that w n = (mod 2n — 3). Note that, modulo 2n — 3, we have 
that 



w r . 



2n-4 

1 - x) XI (-j + = -(2n - 4)! [(* - I) 2 "" 3 

,Y2n-3\ . N ,2n-3/2n-4\ 

where C m := ^y( 2 ^) is the mth Catalan number. The Catalan numbers are integers 
that arise in numerous counting problems. □ 

Remark. It might be interesting to see whether the integers v n j (2n — 3) 1 with i — 1, 2 or 
3 also have a geometric meaning. 

The next result was obtained in collaboration with Alexander Blessing. It estab- 
lisheds parts 2 and 3 of Theorem ^ 
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Lemma 3 For I > we have vi n+ \ = vi n+ 2 (mod 2n) . 

Proof. Since v\ = v 2 — 1 the result is trivially true for Z = and thus we may assume 
/ > 1. We have, modulo 2n, 



2/n-l 



Vin+1 '- 

Furthermore, we have, modulo 2n, 

2ln+l 



J'=0 



2/d 



V/n+2 



Yl (1-j+jx) 

3=0 
lln 

[1 - x)Y[(l - j + jx) 

j=0 
2/n 

[l-x)Y[(l+j-jx) 

3=0 
2/n-l 

vi-z) n (-1 -j+j^) 
3=0 



(1 - x)x JJ(1 -j+jx) 



3=0 



2/d 



;i - - (2ln - 3) + (2ln - j)x) 

3=0 

2/n-l 

;i_. T ) J] (_i)(_i -j+jx) 
3=0 

= v in+1 . 



□ 



This concludes the proof. 

The next lemma generalizes Lemma ^ It implies part 4 of Theorem 
Lemma 4 If k is odd, then vik+(k+3)/2 = (mod (2/ + \) 2 k 2l+2 ). 

Proof. We have v lk+ ( k+3 y 2 = [(1-x) Ylfio 1)fc ((2/ + l)k -j + jx)] x ik+( k +i)/2. The terms in 
the product with j = and j = (21 + l)k lead to a factor of (21 + l) 2 k 2 . The remaining 
terms in the product that are divisible by k lead to a factor k 21 . □ 

3.2 The sequence {v n }^Li modulo primes 

The following lemma will be repeatedly used in this section. 
Lemma 5 Let p be a prime and c an integer. Then, modulo p, 



p-i 

Y\_(^ x — i + c) = 

i=i 



-(x — l) p 1 if p\c, 

-(x + x 2 + . . . + x' 1 '" 1 ) = otherwise. 



Proof. If p\c, then the result is trivial, so assume p\c. We can write 

p— 1 p— 1 



p-i 

i=l 



+ c) = l[t Y[(x-l + c/i) = j^- X 



i=l i=l 



— X 



where we have used that as i runs over 1, 2, • • • ,p — 1, —1+c/i runs over all residues 
modulo p, except for —1. □ 



The second lemma in this section is part 1 of Theorem ^ 
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Lemma 6 For n > 1 we have v n = 1 (mod 2) . 



Proof. Modulo 2, the 2n — 2 terms in the product in Q are alternatingly 1 and x. It 



□ 



thus follows that v n = [(1 — x)x n = 1 (mod 2). This concludes the proof. 

The next lemma together with Lemma El establishes part 5 of Theorem ^ 
Lemma 7 Let p be a prime. Then v p+ \ = v p+2 = 1 (mod p). 

Proo/. Modulo p we have = [(1 - x) ]\^Z\{-1 -j- jx)% P = [(1 - x) (^) 2 ], P , 
where in the derivation of the first congruence we noted that modulo p the jth term in 
the product Q is equal to the (j + p)th and in the second we used lemma |3J Now note 
that 



:i-x) 



X 



X 







n 





X 



X 



E 

k>2 



X 



1. 



xP 



Finally, by Lemma El v p+2 satisfies the same congruence as v p+ i modulo p. □ 
The proof of the next lemma involves congruences for binomial coefficients. In all 
cases these can be found by direct computation, but often it is more convenient to 
invoke a classical result of E. Lucas. Let n > m be natural numbers and write n = 
a Q + a\p + CL2P 2 + ■ — h a s p s and m = bo + b\p + b 2 p 2 + ■ — h b s p s with < aj, bi < p — 1. 
Then Lucas's theorem states that 



(mod p). 



Recall that (?) = if b > a. For example, by direct computation we find that 

;-2) ■ • • (-2p + 1; 



p 2 -2 
2p-2 



1 — p) = — (2p — l)(p — 1) = —1 (mod p) 



([■ 



l---(2p-2) 

means skipping multiples of p). By Lucas's theorem we find that 



p 2 -2 
2p-2 



(p - l)p + p - 2 
1 ■ p + p — 2 



p — 1 
1 



p-2 
p-2 



Likewise we immediately find using Lucas' theorem that, with r 

'2p - L 



•1 (mod p). 
1 and p > 3, 



= 1 (mod p r 



(This identity with r = 2 was proved in 1819 by Charles Babbage. For r = 3 it follows 
from Wolstenholme's theorem, see F.L. Bauer pQ.) 

At various points we use the easy result that ( p_ ) = (— l)- 7 (modp). To see this 
observe that modulo p the entries, except the two outmost ones, in the (p + l)th row 
of Pascal's triangle are zero modulo p. Since each of these entries arises as sum of two 
elements above it in the pth row, the entries in the pth alternate between 1 and -1. 
Similarly one infers that ( P J 2 ) = (— + 1) (mod p). 

For a nice survey of arithmetic properties of binomial coefficients we refer the reader 
to Granville [§]. 
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Lemma 8 Let p be a prime and 2 < I < p + 1. Then v\ p+ i = vi p+2 = —1 (mod p). 

Proof. We have = (1 — x)P l {x) (modp), say, with P l (x) := Yl^lliix — i — l) 21 . 

L J x 1 p 

As before 

1 21 



[1 zi 
(1 - ^-^(l + x + x 2 + . . .) (modp) 

2/ 

[JJ x<-')< (-!)< P) • JJ *« P ] (mod p) 



= x 

L z — ' v ' V i I 

\y — v 

a i(p-l) 
Z-l 

. . . + X ^ a {l-l-i)(p-l) b (i+l)(p-l)-l-l 
i=0 

l-l 

+x i(p ~ 2) ^2 a (i-i-i)(p-i) b (i+i)( P -i)-i + ■■■ ( m ° d p)> 

j=0 

where we have used that 

[21 + n - 1\ 

x'\ 



v fc>0 / n>0 



Thus 



Wip+i = y^g(f-i-i)( p -i)(b(t + i)(p-i)-f - (mod p) 

j=0 

^ /(i + l)(p-l) + Z-2\ 

= 2^a(/-i-i)(p-i) ( 2/ _ 2 l(modp). 



Note that 



2/ \ f -2Z for* = Z-2; 
a ^-^) = ^ _ i _ ij = \ 1 for i = I - 1, 

and, furthermore, 

J (mod p) for i = 0, . . . , p — I — 2; 
a(,_i_ iK p_i) = | _ 1 ^ mod ^ f OT i = p - Z - 1. 

The remainder of the proof requires a separate discussion for I < (p — l)/2, for I = 
(p + 1) /2, (p + l)/2 < I < p, for / = p and for I — p + 1. 
The case Z < (p — l)/2 : 
Note that 

(i + l)(p-l) + Z-2\ [(i + l)p-i-Z]---[(i + l)p-i + Z-3] 



2/ - 2 J 1 • • • {21 - 2) 

(mod p) for i = 0, . . . , Z — 3; 

( ~1?;.^7_ 2 2| 2) = 1 (mod p) for 2 = I - 2; 

( " 2 ?;.^7_ 2 2| 1) = 2Z - 1 (mod p) fori = Z-l. 
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Hence in the sum for only the terms i = I — 2 and i — I — 1 contribute, that is: 

vi p+ i = -2M + 1 • (21 - 1) = -1 (mod p). 
The case Z = (p + l)/2 : 
Here we find that 



(z + l)(p- 1) + Z 
2/ - 2 



1 (mod p) for i = I — 2; 

(mod p) for < i < I - 1, i ^ I - 2. 



Hence in the sum for Vi p+ \, only the term i 
Vi p+ i = —21 • 1 = — 1 (mod p). 
The case (p + l)/2 < Z < p : 
Here we compute that 



Z — 2 contributes and we infer that 



(i + l)(p-l) + Z 
2/ - 2 



(mod p) for z = p — I, . . . , I — 3; 

—I (mod p) for i = p — I — 1; 

2 — Z (mod p) for i = I — 2; 

(1 - 2/)(/ - 1) (mod p) fori = Z-l. 



Hence in the sum for v ; p +i, only the terms i = p — / — 1, z = / — 2 and i = l — l contribute. 
That is: v ip+1 = -1 ■ (-/) - 21 ■ (2 - /) + 1 ■ (1 - 21) (I - 1) = -1 (mod p). 
The case I — p : 
Note that 



a 



(p-i-i)(p-i) 



(-1) 



p— 1— i 



2p 
p — 1 - 



(mod p) for z = 0, . . . , p — 2; 

1 (mod p) for i = p — 1, 



while for z = p— 1, ( ( ' i+1 ^2p- 1 2 +P 2 ) boils down to (2 P _2) = — ^ ( m °d p) ( see the preamble 
to this lemma). Hence = — 1 (mod p). 
The case Z = p + 1 : 
Here we find that 



a 



(p-*)(p-i) 



(-1) 



2p 
p- 



-2 (mod p) for z = 0; 

(mod p) for z = 1, . . . , Z — 4; 

1 (mod p) for z = Z — 3; 
-2 (mod p) for z = Z — 2; 

1 (mod p) for z = Z — 1. 



and 



{0 (mod p) for i = 0; 

1 (mod p) for z = Z - 3; 

1 (mod p) for z = Z — 2; 

(mod p) for z = Z — 1. 

Hence v lp+1 = —2-0 + 1-1 — 2-1 + 1- = — 1 (mod p). 

This proves that in all cases, v lp+ i = — 1 (mod p). By lemma El the same is true for 

Vlp+2- D 

The cases Z = p and Z = p + 1 in the latter proof can be proved more succinctly as is 
done in the proofs of Lemma |Hl respectively Lemma ITU1 



Lemma 9 Let p be a prime. Then v p 2 +1 = v p 2 +2 = —1 (mod p). 
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Proof. Note that, modulo p, the integer v p 2 + % is congruent to 



p-i 

(1-x) "[[(-!- j-jx] 
3=1 



2p 



xP 



p-1 

3=1 



P n2 



p-1 



On proceding as in the previous proof we find that 

2 / \ 2 



V p 2 + l 



y p -y 
i-y 



J yp 



i-y 



=-l(modp). 



fc>l 



Finally, by Lemma El "y p 2 +2 satisfies the same congruence as v p 2 +1 modulo p. 
Lemma 10 Let p be a prime. Then v p 2 +p+1 = v p 2 +p+2 = — 1 (mod p). 
Proof. We have 



□ 



v p 2 +p+i 



(1 -x)(x + x 2 H h/ ) 



,p-l\2p+2 



(mod p) 



1 - x)(x + x 1 + • ■ • + x v ~ x f{x v + • • • + X 



rP(P~ !h 2 



(mod p) 



(x - x p ) (x + x 2 H h x p ' x ) ( ^2 xkV ) ] 2+ ( mod P) 

k=l XP P 



- [< 



a; 2 + • ■ ■ + x p - x p+1 



oo 

[$> + !)* 



(fc+i)p 



- x 2 ^ 1 ) ^(ifc + l)x^ 

fc=0 

-1 (mod p). 



fc=0 



Finally, by Lemma El v p 2 +p+2 satisfies the same congruence as v p 2 +p+1 modulo p. □ 

The next lemma establishes part 6 of Theorem ^ 
Lemma 11 If p is an odd prime, then Vp+3 +i = (mod p) for i — 0, . . . , (p — 3)/2. 

Proof. In case i — 0, the result follows by Lemma Q] so assume that z > 1. On using 
that modulo p the jth term equals the (j + p)th term, we find that, modulo p, 



V p+3 



+i 



p-1 2i 

4i 2 (l - x) Y[(2i — j + jx) Y[(2i — j + jx) 

j=l i=i 



. x (p+l)/2+j 



On invoking Lemma El and noting that p > + i we infer that 



2/ 



2; 



4i 2 (x p - x) \\{2i - j + jx) 

i=i 



c (p+l)/2+i 



4* 2 JJ(2i- j+ jx) 

3=1 



Since deg^ rj 2 ^ 1 (2i — j + jx)j = 2i and 2i < v -^- + i, the result follows. □ 
The next lemma will be used in the proof of Lemma [E 
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Lemma 12 Define A r (x) and B r (x) recursively by 



A (x) = 0, A r+1 (x) = ( x + ... + x p ~ l Y - A r (x); 
B (x) = 0, B r+1 (x) = -( x + ... + x p ~ l ) r - B r (x). 

Put f r (x) = (x - 1)(1 + x p + . . . + x p (p-V)(x + ... + x p ~ 1 ) r . Then 

f r ( x ) = (-l) r (x - 1)(1 + x p + . . . + x p{p - 1] ) + x p2 A r {x) + B r {x), 
where, for r > 1, the degree of B r (x) equals (r — l)(p — 1). 
Proof. Easily follows on noting that 

f r+1 (x) = (l + ... + x p - 1 )/ r (x)-/ r (x) 

= (X p2 -1)(X + ... + X p - 1 ) r - f r (x). 

The next lemma is part 7 of Theorem ^ 

Lemma 13 Letp be an odd prime. Suppose thatO <i< (p— 3)/2 and (p— 1)/2 < I < p. 
Then vi p+{p+3)/2+i = (mod p). 

Proof. Write / = (p — l)/2 + k (thus k assumes the values 0, . . . , (p + l)/2). Put 
P2i{x) = (2i) 2l+2 Ylf=i(^ ~ J ' + 3 X )- Note that the degree of this polynomial equals 2i. 
Proceding as in Lemma El we infer that, modulo p, 

\ X -l)( X+ ... + X P-Y +1 P2i(x) 



x lp+(p+l)/2+i 

x - \){x p + x 2p . . . + s p(p - 1) )(a; + . . . + a* -1 ) 2 *^^ 



t .lp+(p+l)/2+i 



(6) 



We consider the case I = p (that is k = (p + l)/2) first. Then 

(z - l)(x p + x 2p + . . . + x p(p - 1} ) 2 (x + . . . + x p - 1 )P 2i (x) 
(x p - x)(x p + x 2p + . . . + x p{p - 1) ) 2 P 2i {x) 



c lp+(p+l)/2+i 



_ x lp+(p+l)/2+i 



where we used the observation that the polynomial in brackets has the form 

P2i(x)(c kp x kp + c kp+1 x kp+l ), 

k 

and that 1 + 2i < (p + l)/2 + % < p - 1. 

Thus we may assume that I < p — 1. Notice that /p+(p+l)/2 + z < p 2 . Furthermore, 
by Lemma IT2l we have [B 2k (x)P 2i (x)] x {i-i)p+( P +i)/2+i = 0, as 

2k - 1 + 2i < (I - l)p + V -^- + %. 

Using this, (jOJ) and Lemma IT2l we infer that 

(x-l)(x p + ... + x p(p - iy )(x + ... + x p - 1 ) 2k P 2i {x) 
{x - l){x p + ... + x p2 ){x + ... + x p - 1 ) 2k P 2i (x 
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-ip+(p+l)/2+i 



_ x lp+(p+l)/2+i 



(x- 1)(1 + . 
(x-l)(l + . 
(x-l)(l + . 

This concludes the proof. 



. +x P(P-l))( x+ ... +x P-l)2fcp 2i ( x ) 

.+x p ^)P 2t (x)+B 2k (x)P 2l (x 

. + x p(p-l))p 2i ( x ) 



x (l-l)p+(p+l)/2+i 



c (l-l)p+(p+l)/2+i 



_ x (l-l)p+(p+l)/2+i 

= 0. 



□ 



A further question concerning the distribution of v n modulo primes is how frequently 
certain residues appear. For example, is it true that the zero have density 1 ? Is it true 
that the non-zero entries are equidistributed ? Questions like this can be answered for 
the middle binomial coefficient ( 2 fc fc ), see e.g. (211111 El- The following lemma suggests 
that perhaps techniques from the latter papers can be used to investigate this issue. 

Lemma 14 We have Vi + ^k = v 2+ zk = k+iCk) (mod 3) and v 3 k = (mod 3). 

Proof. By Lemma Q we have = (mod 3). By Lemma El we have t>i+3fc = 
v 2+ 3k (mod 3). We have, modulo 3, 



V 2+ 3k 



[1 - x)(l + x 
This concludes the proof. 



l+6fc 

3=0 

2k 



r l+3k 



(1 -x)(-x(l +x)) 2k x 



x l+3k 



2k 
k 



2k 
k-1 



1 



k + l\ k 



2k 



□ 



3.3 The sequence {v n }™ =1 modulo prime powers 

The proof of the next lemma was kindly communicated to us by Carl Pomerance. 

Lemma 15 The polynomial Y%^~{ix —i + j) (mod p l ), as a polynomial in x, depends 
only on the class j (mod p) (ie. replacing j by j + kp would yield the same result). 

Proof. Let fj(x) = Yli=o + j)- ^ P\h then there are p r ~ l factors divisible by 
p and p r ~ 1 > r, so that fj(x) = (modp r ). So assume p \ j. Let k be the in- 
verse of j, so jk = 1 (modp r ). Then modulo p r , we have fj(x) = j pT fi(x) (since 
the expression ik runs over a complete residue system modulo p r as i runs). Now say 
j = ji (mod p) , say ji = j + kp. Using induction with respect to r one then easily sees 
that j p = (j + kp) pr = j pr (mod p r ), and we are done. □ 

Proof of Theorem^ Part 1. We have 

p-i 



V P+3 
2 



P 



P 



X 



X 



\\(P- 3+ 3%) 



p-1 



3=1 
p-1 



p—1 p— 1 

Y[(-j +jx) + (i- x)p Y[(-j + j x ) 



3=1 



k=l i=i 



p-1 



+ (l-x)p 2 Y Y[(-j+3 x ) p-i (mod p 5 ) 



l<fc<r<p-l J =1 

j¥=k,r 
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p-1 



= -p 2 ( P -i)\ (x-iy + ix-iy-^y^-r 

+ {x-l) p ~ 2 p 2 Y p-i (mod p ; 

l<fc<r<p— 1 

ee -p 2 {(p -l)\}[(x- If] ^ (mod p 5 ), 
= ~P 2 {(P - I)'} (A) (-1)^ (mod p 5 ), 
ee 2p 3 {(p -1)!}(1 -p) 



p — 1 

p-i 

2 



1) 2 (mod p ), 



where we used that ^fc=i 1/^ = (mod p 2 ) (this is Wolstenholme's theorem ^21 Theo- 
rem 115]) and Yli<k<r<p-i ^/kr = (mod p). To see the latter congruence note that 



l<fc<r<p-l j=l 



j.p-3 



rP-1 



j.p-3 



(mod p). 



Now it is an easy consequence of Eisenstein's congruence (1859), see O Theorem 132], 
which states that 

2P- 1 -1,11 1 . J , 

= 1 + - + - + •■■ + (mod p), 

p 3 5 p — 2 

that (see ibid. Theorem 133) (( p Jn/ 2 )( — l) 2 ^ - = 4 P_1 (modp 2 ). (Indeed, by Morley's 
congruence (1895), cf. j3j, this congruence is even valid modulo p 3 .) We thus finally 
infer that v p+3 ee 2p 3 (l — p){(p — 1)!}4 P_1 (mod p 5 ), which of course implies that Vp+a = 
—2p 3 (mod p 4 ). 

Part 2. We have the formal series identity 



(1-x) 



2r 



k=0 



k + 2r - 1 
2r - 1 



x 



Note that 



Using the latter congruence we find that, modulo p 2r + 3 



(1-23)2,-! 

= -s — n rr (mod »). 

22r-i(2r- 1)! v ^ ; 



17p+3 



+rp 



(2r+l)p 

(1-ar) JJ ((2r + l)p-j+jV 

3=0 

2r+l 



P 



2r+2 



p+i_ 



(2r+l)p 



(l-x)Y[(2r + l-j+jx) JJ (-j + jx) 



j=0 



2r 

(2r + 1) 2 p 2r+2 [( hrx j ) {x - l) (2r 



3=0 
P\j 

+l)p-2r 



P+i 



3=0 
2r 



a ,(p-l)/2+rp 



(2r + l)V^ 2 \(Y,h 



X J )(x 



^(2r+l) P y^ 



3=0 



k=0 



k + 2r-l 
2r - 1 



x 



x (p-l)/2+rp 
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c (p-l)/2+rp 



- -(* + 1) V« [( £ 6, r ^) J] x'-(-l)' p (* + * " ') ; 

1=0 ^ ' j=0 ^ ' 

( 2r + l)2 p 2r+2 r . , 2r 

IT-M Z (- 1 )'EM(l-2j)) 2 ,-i 

' ;— n V / i— n 



2 2 '- 1 (2r -« , 

_ (-l)^(2r + l)V 2A * ' 
" 2 2 -i(2r-l)! J)V 

where in the one but last step we used the identity 



which is obtained by comparing the coefficient of x r of both sides of the identity 

(1 — — x) 2r+1 = (1 — x) 2r . This finishes the proof. □ 

3.4 The sequence {v n }™ =1 modulo powers of two 

Before we can consider the sequence modulo powers of two we need some preparatory 
lemmas. 

Lemma 16 If j is odd, then Yli=o {i x ~~ i + j) 2 = x<iq (mod 2 q ). 

Proof. By induction with respect to q. For q — 1 the result is obvious. Assume the 
result is established for 1 < q < q 1 . We write 

2<n+ 1 -l 2^1-1 2«i +1 -l 

Yl (ix-t+j) 2 = Y[(tx-i + j) 2 Yl (ix-i + j) 2 = P 1 (x)P 2 (x) 1 

i=0 i=0 i=2H 

say. Note that Pi(x) = P 2 (x) (mod 2 qi ). The induction hypothesis thus implies that 
we can write P\(x) = x 291 + 2 qi fi(x) and P 2 (x) = x 2n + 2 qi f 2 (x). Since (ix — i + j) 2 = 
((z + 2 qi )x - (i + 2 qi ) +j) 2 (mod 2 9l+1 ), it even follows that P x (x) = P 2 (x) (mod 2 qi+1 ), 
from which we infer that f\(x) = f 2 (x) (mod 2) and hence f\(x) + f 2 (x) = (mod 2). 
It follows that modulo 2 qi+1 the product under consideration equals 

P l (x)P 2 (x) = (x 2n +2 qi f 1 (x))(x 2 " 1 +2 qi f 2 (x)) =x 2n+1 (mod 2 9l+1 ). 

This concludes the proof. □ 

In the course of the above proof we have showed that 

29-1 2« +1 -l 

Y[{ix-i+j) 2 = Yl (ix-i + j) 2 (mod2 9+1 ). 

i=0 i=2« 

The next result shows that the same identity holds true for the 'square roots'. Using 
this the 'square root' of the left hand side of Lemma Hoi can be computed (Lemma I18|). 
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Lemma 17 Let j be odd and q>2. Then 

21-1 2i +1 -l 

Y[ (ix-i+ j) = ] [ (/./• - / 4- j) (mod 2 q+1 ). 

i=0 i=2i 

Proof. It is an easy observation that, modulo 2, we have for < k < 2 q — 1 that 

2«-l 



jr { x it k is even. 



a=0, a^k 

Using this identity we find that, modulo 2 q+1 , 

29+1-1 2«-l 



[ [ (ix-i + j) = ]J(ix-i+j + 2 q (x - 1)) 

i=0 

21-1 

l[(ix-i + j) + 2' ! (x-l)J2 II fa-i+j) 

i=0 k=0 i=0, i^k 

21-1 21-2 21-1 

H(ix-i+ 3 ) + 2%x - 1) (x 2 "- 1 i + ^ 9 " i_i 1 

i=0 2\k 2\k 

21-1 

\\ (ix-i + j) + 2 q [x - l){x 2 "~ 1 2 q - 1 + x 2 "' 1 -^- 1 ) 

i=0 
21-1 



i=29 i=0 

21-1 21-1 21-1 



i=0 

This finishes the proof. □ 
Lemma 18 Let j be odd and q>3. We have 

21-1 

(mod 2 q ). 



Y[{ix-i + j) = x 2 " 1-2 2 q ~ 1 (x A + x 3 + x + 1) + x 2 



i=0 

Proof. Similar to that of Lemma Hoi but with the difference that instead of the equality 
Pi(x) = Pzix) (mod 2 q+1 ), we use Lemma IT71 □ 

Remark. By Lemma El it suffices to work in the proofs of Lemma El El an d El with 
j — 1- 

The next result established part of parts 8 and 9 of Theorem ^ 
Lemma 19 For q > 1 we have V21 = — 1 (mod 2 q ). 

Proof. Put P q (x) = (1 — x) rijlo 3 ( — 3 — j + jx). We want to compute the coefficient 
of x 2<i ~ x in P q (x) modulo 2 q . On invoking Lemma IT^l one finds that 

P q (x)(l + 2x){2 + x) = x 2 \l - x) (mod 2 q ), 
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from which we infer that 

oo 

P q (x) = x 2q ' q {l - x) J^(-2) V ^(-2)' ? - 1 - r x r (mod 2 q 



q-l 



k=0 r=0 

2q-2 2,-1 

X 2 "- q (l ~X)J2 a mX m = X 29 - Q b mX m (mod 2"), 

m=0 m=0 



where 



J -(-2)' ? - 1 - m /3 ifO<m<g-l; = / -(-2)^ 1 - ifO<m<g-l; 

\ -(-2)- q+1+m /3 if q < m < 2q - 2 m ~ \ (-2) m -<? 



Thus w 29 = 6,_i = -1 (mod 2 q ). 



if g < m < 2q — 1. 

□ 



Recall that we defined vq = — 1. The reason for this is that this definition allows us 
to also formulate the next lemma, which together with Lemma I2S1 is part 8 of Theorem 
□ with j = 0. 

Lemma 20 (Periodicity.) Suppose that i, k > 0. We have v^i+i = Vi (mod 2 q ). 
Proof. First assume that i > 2. Note that 

2«-l 2i-3 



Vk2i+i 



[I - x) H (2i - 3 - j + jx) 2k l[ (2i - 3 - j + jx) 

j=0 j=0 



r k21 + i-l 



fmod 2 q ). 



By lemma ITol the first product equals x k2Q mod 2 q . Thus 

2i-3 

v k 2i+i= (1 - x) JJ(2i - 3 - j + jx) 

j=0 



u,- (mod 2 9 ). 



In order to deal with the case i — 1, we note that, using Lemma El tto+i = ^29+2 = 
t> 2 = t>i (mod 2 q ). In case z = one finds proceding as above that, for k > 1, = 
t>29 (mod 2 9 ). On invoking Lemma ITU1 it then follows that v^t = V21 = Vq (mod 2 q ). □ 

The next result yields a part of part 9 of Theorem ^ 

Lemma 21 Suppose that q > 1. T/ien t> 29 -i = 2 9_1 — 1 (mod 2 9 ). 

Proof. Similar to that of Lemma For q < 3 one verifies the claim numerically. So 
assume g > 4. We want to compute the coefficient of x 2Q _1 in P q _i(x) modulo 2 q . On 
invoking Lemma ITH1 one finds that 



P q -i(x)(l + 2x)(2 + x)= x 29 '- 2 2 q -\x A + x :i + x + 1) + x 2 (mod 2 q ) 



whence P q -i(x) = x 2 "' 1 -"' 2 ( p-^x 4 + x 3 + x + 1) + x 2 J (mod 2 9 ). (Note 

that the assumption q > 4 implies that 2 9_1 — q — 2 > 0.) Thus modulo 2 9 the coefficient 
of x 2q _1 , that is U29-1, equals 



U 2 «-i = + ^-S + &,-2 + &, + &,+l) 
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= -1 + 2 



9-1/ 



-4 + 2-2 + 1} 



2 q - 1 - 1. 



This completes the proof. 



□ 



The next result with % = 0, 1 and 2 yields a part of part 9 of Theorem It also 
yields part 10 of Theorem ^ 

Lemma 22 For i > and q > 2 we have v 2 i-i+i = u j + 2 9-1 (mod 2 9 ). 

Proof. For g = 2 one checks the result numerically and so we may assume q > 3. For z = 
the result follows by Lemmal^Tl Note that, a priori, u 2 «-i+j = t>j (mod 2 9 " 1 ) and so either 
f 2 <j-i+j = Vi (mod 2 9 ) or u 2 «-i+j = Vi + 2 q ~ 1 (mod 2 9 ). Let us first assume that i > 2. The 
idea of the proof is to use Lemma IT7I to write v 2 q-i+i = v i + 2 q ~ 1 [f q> i(x)] x2 q-i + i- 1 (mod 2 q ). 
An easy computation then shows that [fq,i{x)\ w-i+i-i is odd, thus finishing the proof. 
More precisely, one first notes that, modulo 2 q , 



21-1 



v 2i- 1 +i 



29+2i-3 



[l-x) ]J(2i-3-j+jx) ]"[ (2i-3 



x 21- 1 +i-l 



IX 



3 + 3%) 

3=0 i=2« 
2«-l 2i-3 

1-x) J](2z-3-j+jx) J](2z-3-j+j: 

j=o i=0 

2i-3 

1 - x)x 2q ' 1 - 2 {2 q - 1 {x A + x 3 + x + 1) + x 2 ) Y[ (2i - 3 - j + jx) 

2i-3 

^i + 2^ 1 (1 -x)(x 4 + x 3 + x + l) JJ(2i - 3 - j + jx) 

j=0 



-i 



^1 — x)(x 4 + x 3 + x + l)x 



i-1 



Since [(1 — x)(x 4 + x 3 + x + l)^ -1 ]^! = [(1 — x)(x 4 + x 3 + x + l)] x .2 = —1 is odd, the 
result follows for j > 2. On combining this result for i — 2 and Lemma E3 we find that 
U29-1+1 = v 29-1+2 = v 2 = v i (mod 2 9 ). Thus the result follows for every j > 0. □ 

Using induction and the latter lemma one then easily infers the following result 
which, together with Lemma EH1 gives part 8 of Theorem ^ 

Lemma 23 (Equidistribution.) Let q > 1. For even/ odd integer a there are precisely 
two integers 1 < ji < J2 < 2 17 such that Vj 1 = a (mod 2 q ) and Vj, 2 = a (mod 2 9 ). 



3.5 On a result of Paolo Dominici 

Let Sk(xi, . . . , x r ) denote the fcth elementary symmetric function in r variables, i.e. 
Si(xi, . . . , x r ) — x\ + . . . + x r , S2(xi, . . . , x r ) = X\X2 + X1X3 + . . . + x r _ix r , etc.. Paolo 
Dominici [7j states the following result for v n without reference. 

Theorem 3 For 1 < i < 2n — 4 we put Ui — i/ (2n — 3 — i) . Then 

V n = (2n - 3) 2 (2n - 4)!{S , „,_ 2 (2/l, • • -,y2n-i) ~ Sn-lG/l, • • • ,V2n-i)}- 

We will now derive this result from (JHJ). We need two lemmas 
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Lemma 24 Let Li(x), . . . , L r (x) be linear polynomials, then 



1 d m 
ml dx m 



{Li(x) ■ ■■L r {x)} = S ri 



|M,...,|M |ilW ... Ml , 



Another observation we need is the following: 

Lemma 25 Let xi, . . . , x r be distinct non-zero elements such that x\ • x 2 • • • x r — 1 and 
{xi, ...,x r } = {i, . . . , j-}. Then S r ^ k {x u ...,x r ) = Sfc(xi, . . .,x r ), with 1 < r < k. 

Proof. Note that S r _k(%i, ■ ■ ■ , %r) — ^k(^~, ■ ■ ■ > x ) x i ' ' ' x r — Sk(%i, ■ ■ ■ , x r ), where in 
the derivation of the first equality we used the assumption that ^ and in that of 
the second the remaining assumptions. □ 

Corollary 1 For 1 < k < 2n - 5 we have S k (y l7 . . . , y 2n -±) = S 2n -A-k{y\, V2n-d- 

Proof of Theorem^ Put P n (x) = rij=i~ 4 (2 n ~~ 3 — j + jx). By definition we have 

2n-3 



[l-x) H(2n-3-j+ji 

3=0 



(2n-3) 2 (l-x)P n (x) 



Thus, 



v n = (2n - 3) 2 {[P n (x)} xn -2 - [P n (x)] x n-a}. 



(7) 



On noting that [P n (^)]a; 



J__£ 

ml dx 



^P n (x)\ x=0 , we obtain on invoking Lemma EH that 



[Pn(x)] xm = (2n - 4)\S m (yi, . . . , y t , . . . , y 2n -i)- 
Combining (JSJ) with (JJJ) yields that 

w n = (2n - 3) 2 (2n - 4)!{S , n _ 2 (2/i, . . . , y 2 n-4) - S n - 3 (yi, • • • , y2n-4)}, 
or, on invoking Corollary Q 

v n = (2n - 3) 2 (2n - 4)!{S , n _ 2 (l/i, • • - ,yan-4) - Sn-ifai, • • • ,3/2n-4)}, 
This concludes the proof. 

From the above proof we infer that we may alternatively define v n by 

2n-3 

JJ(2n-3-i+jx) 

j=0 



(8) 



□ 



(9) 



We leave it to the reader to use the observation that P(x) := n 7 -=i ~ 3 — j + jx) is 
selfreciprocal, i.e. satisfies P{l/x)x 



2m-4 



P(x) to infer (jUJ) directly from (j3*jl. 



Many of the congruences can be also proved using Theorem El As an example we 
will show that if p is an odd prime, then v 3 ( p+ iy 2 = — 81p 4 (mod p 5 ). This is the case 
r = 1 of part 2 of Theorem |21 
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Proof of part 2 of Theorem\^in case r — 1. Set n = 3(p+l)/2. Note that (2n — 3) 2 (2n — 
4)! = — 18p 4 (modp 5 ). It thus remains to be proven that the expression in braces in 
Theorem El equals 9/2 modulo p. It turns out to be a little easier to work with Wi = —y^. 
Note that (— l) r S r (wx, . . . , ii^p-i) = S r (yi, . . . , jfep-i). We have = — 3p) for 
1 < i < 3p— 1. Thus W7 P = —1/2, u^p = —2 and the remaining Wi satisfy Wi = 1 (mod p). 
Hence S r (wi, . . . , ^3 P _i) = SV( — 1/2, —2, 1, 1, . . . , 1) (mod p), where 2 < r < 3p — 1. In 
the symmetric function S r (zi, . . . , z 3p _i) there are ( 3p ~ 3 ) terms containing neither Z\ nor 
z 2 . There are ( 3 ^Zf) terms containing z 1; but not Finally there are ( 3 ^Z 2 3 ) t erms 
containing both z\ and z 2 . It follows that, modulo p, 

(-!/*(», . . . ,^0 , -2, 1, .... i) = ( 3p ; 3 ) - (2+ i) ( 3 r " :i 3 ) + ( 3 ;; 3 ) . 

Modulo p we have 



•l) n {5 , n _ 2 (?/l, • • • , V2n-A) ~ S n -l(yi, 2/2n-4)} 

'3p -2\ 5 /3p - 2\ /3p - 2 
n-l J ~2\n-2 J + \n - 3 

P-2 \ t P-2 \J P -2 
\n — p — 2 J \n — p — 3 
p-2 



n — p — 1 



21 P ~ 2 J+4, 

n — p — \ J \n — p — 2 



p-2 
n — p — 3 



p 

n — p — 1 
p-2 
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p-2 
n — p — 2 



-9. 

— p — 2 
This completes the proof. 



■l)»9(n-p- !) = (-!)»-. 



p-2 
n — p — 2 



□ 



4 Asymptotics 

Given a sequence of coefficients, there are many things we would like to know about 
it. Apart from the search for a generating function and for a recursion formula, an 
interesting question is the asymptotic behaviour. We remind the reader that candidate 
Fourier series for modular forms of weight 2k for SL(2, Z) must have coefficients growing 
like n 2k ~ l (and n k for cusp forms). 

In our case, without prior knowledge of the alternative definition (jSJ), we only man- 
aged to compute the first 80 values of v n using the Schubert package for intersection 
theory ^Hj or the first 225 values using the rational function (with the dummy variables 
Wi). Numerically, it is readily seen that the leading term for the v n is e 2nlogn . As this 
is strongly reminiscent of the behaviour of (2n)! = exp(2nlog(2n) — 2n + |log2n + 
|log27r + 0(~)), we rather study the behaviour of log and find now the leading 
term to be 2n. Subtracting it, we find the next-to-leading term to be — 41ogn, easily 
verified by applying nd n (ie. taking subsequent differences and multiplying by n). The 
next term is a constant, C = —5.62..., which we find difficult to recognize. We have 
learnt from Don Zagier a smart technique which enables to determine a large number of 
digits of C; we present it below under the name of asymp^ trick. 
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4.1 The asympfc trick 

Assume we are given numerically a few hundred terms of a sequence s = {s n } ne N which 
we believe has an asymptotic expansion goes in inverse powers of n, ie. 

c i c 2 

Sn ~ C H h — + . . . . 

n rr 

Goal : determine the coefficients q numerically. 

Trick : Choose some moderate value of k (say k — 8) and define a new sequence as 
—A k N k s, where A is the difference operator (Au) n = u n —u n -\ and N the multiplication 
operator (Nu) n = nu n , i.e., 

(fc) ( "0^ / -\k 

For n large we have (assuming the above asymptotic expansion for s itself) 



U k c fc+i / n fc (( fc + x ) c fc+2 - (T)cfc+i) 



4 fc) = c + (-ir^+(-ir^ — ' ^-^^ 

Thus, while s n approximates Co only to within an accuracy 0(n _1 ), approximates 
it to the much better accuracy 0{n~ k ), so we obtain a very good approximation for cq. 
Call this operation asymp^. The further coefficients q are then obtained inductively: if 
Co, ... , Q-i are known to high precision, we get q by applying asymp^ to the sequence 
n l {s n - c Q-i/n^ 1 ) = d + Ci+i/n H . 

The crucial point in the success of asymp/; is that the operator A fc sends n k to k\ 
and kills polynomials of degree < k, so that all the intermediate terms of the expansion 
of s n between c and c^n~ k disappear. 

Variants of asymp^ allow one to deal for example with asymptotic expansions of the 
form 

(I) s n ~ Alogn + c + ci/n + c 2 /n 2 H 

(II) s n ~ Bn + Alogn + c + cx/n + c 2 /n 2 -\ 

(III) s n ~ An x (l + ci/n + c 2 /n 2 + • • •) 

In case (I) we can apply asymp fc to the sequence n(s n+ i — s n ), which has the form 
A + c[/n + c' 2 /n 2 + • • •, to obtain A to high precision, after which we apply the original 
method to {s n — Alogn}. In case (II) we apply asymp^ twice to As to get B and A, 
and then subtract (our approximation for) B log n + A from s n and apply the standard 
version. For case (III) we can either look at {logs„} and apply variant (I) or else apply 
asymp fe to {n(s n+1 /s n — 1)} to get A and then apply the standard method to {s n /n x }. 

Remark 1. Applying the operation asymp fe with suitably chosen k gives a rapidly con- 
vergent sequence . To estimate how many decimals are probably correct, we look 
at some relatively widely spaced elements of this sequence (e.g., the terms si^ with 
n = 300, 400, 500 if we know 500 terms of the sequence s) and see how many of their 
digits agree. 
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Remark 2. One also has to experiment to find the optimal choice of k. Typically one 
uses k = 5 if one knows 200 terms of s and k — 8 if one knows 1000 terms. This suggests 
that perhaps k « log N is a good choice for a generic sequence with N computed terms. 

Remark 3. The asymp^ trick was first described in a paper of Zagier [21]. Here he con- 
siders the Stoimenov numbers £r> which bound the number V(D) of linearly indepedent 
Vassiliev invariants of degree D. Stoimenov himself thought that £d behaves 'something 
like D\/1.5 D \ Calculating the values up to D = 200 and applying a variation of asymp^ 
suggested an asymptotic formula of the form 

with C ~ 2.704332490062429595, C x « -1.52707 and C 2 « -0.269009. Subsequently 
Zagier was able to prove this with explicitly computable constants Cj. In particular, 
C = ^v^Tr^/V 2 / 12 , which agrees to the accuracy given above with the empirically 
obtained value. 

4.2 Application to the asymptotics of v n 

In our case of sequence v n of lines in a hypersurface of P n , the coefficients c =: C is 
difficult to recognize, but all other coefficients, c\, c%, . . . are rational numbers which we 
easily recognize from a sufficient number of digits. Once the first few rational coeffi- 
cients have been found and the corresponding terms subtracted from the sequence s, 
the constant term C can be obtained with 30 digits, say. This is enough to feed to 
the PARI software and apply the function lindep( [C, 1 , log(Pi) ,log(2) ,log(3)]) to 
find a rational linear combination of C in terms of a given basis (educated guess). The 
result, equivalent to (J3J is: 

log (^ = 2 "- 41og " + c + ^ + i^ + -- (10) 

where C := — 3 — log7r — | log |. In the appendix we present a proof by Don Zagier of 
this asymptotic formula. 

5 Comparison with two other sequences 

As a matter of curiosity, we now compare our results so far with similar results from 
two other sequences of enumerative geometry. We shall see that the first case has quite 
similar features, while the second case is more intricate. 

5.1 Numbers of plane rational curves 

One sequence of integers from enumerative geometry is na, the number of plane rational 
curves of degree d through 3d— 1 points in P 2 . Kontsevich's recursion formula ^1] reads 



n d = 

k=i 



Til = 1. 
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The result is n\ — 1, n% — 1, n% — 12, etc, ie. there is 1 line through 2 points of the 
plane, 1 conic through 5 points of the plane, 12 cubics through 8 points of the plane, 
etc. 

We can similarly draw tables of rid mod k for any integer k. The results (in the same 
convention as before) are: 

* k = 2: both rows vanish (except first two values), ie. all rid are even. 

* k = 2 l : all rows are 0, ie rid = mod 2 l for n > I + 1. 

* k = 3: n 3 d = mod 3, n 3 d+2 = 1 mod 3, n 3 d+i = alternating 1 or 2 mod 3 because 
n M+2 = 4 mod 6. 

* /c = 5: n d = mod 5, for d > 8. Idem for k = 25 (d > 23) 

Because of this big symmetry for low primes, most non-primes will yield constant or 
regular rows (ie rows repeating when shifting horizontally). The only non-obvious case 
is k = 26, where there is a shift by 8 (because k = 13 shifts by 16) and rows 4,6 alternate 
with 0. 

Further, we only found three primes with regular features: 

* k = 7: all rows are regular (repeat when shifted horizontally by 4), rows 5 and 7 
are 0. 

* k = 13: idem, shift by 16, no row. 

* k = 19: idem, shift by 12, no row. 

* k = 5, 11, 17, 23, 29: these primes give almost-0 rows (ie. rid = except for a finite 
number of d). 

We have not attempted to prove these observations. 
5.1.1 Asymptotics 

We now turn to the asymptotics of the sequence rid for d — ► oo. Di Francesco and 
Itzykson proved [Ul Proposition 3] that 



as d tends to infinity, and found the approximate values A pa 0.138 and B pa 6.1 for the 
constants A and B. Assuming a full asymptotic expansion 



and applying variant (II) of the asymp fc trick to \og(rid/{3d — 1)!), we obtain the much 
more accurate approximations 

A pa 0.138009346634518656829562628891755541716014121072, 
B pa 6.0358078488159024106383768720948935 , 

as well as the further values B 1 pa -2.2352424409362074, B 2 pa 0.054313787925. Unfor- 
tunately, we are not able to recognize any of these apparently irrational numbers, e.g., 
PARI does not see in log A and logi?o a linear combination of simple numbers like 1, 
log 2, log 3, log7r, Ti and n 2 . 
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5.2 Numbers of rational curves on the quintic threefold 

The other sequence we now introduce for the purpose of comparison is qa, that of holo- 
morphic rational curves of degree d embedded in the quintic Calabi-Yau threefold. These 
are the 'instanton numbers' of Candelas et al |S]. They are defined by the following line: 

5 + Vg„n 3 -^=f^^ -i— -_ = 5 + 2875 g + 4876875 g 2 + . . . , (11) 

1 — q n \x dq J (1 — 5 b x) y {xy 

where q(x) = x e^ 1 ^ = x + 770 x 2 + . . . is the "mirror map" and its inverse is x(q) = 
q — 770 q 2 + . . .. The functions yo and y\ are solutions of a Picard-Fuchs differential 
equation and are given by 

, \ v-(5n)! n . „ . , f(5n)\ ^ A n 

w>(*) : =2^^r x and = Lhr 5 L- x 

n>0 n>0 \ j=n+l J J 

When computing the numbers q^, the longest step is without doubt the inversion of the 
series q(x). The first few values are q\ = 2875, q 2 = 609250, etc. 

We can again draw tables of q^ mod k for any integer k. The main results (in the 
same convention as before) are: 

* k = 2: second row is 0, ie. q 2 d are even. 

* k = 4, 8, 16: last row is 0, ie. q 2 i d = mod 2 l {I < 4) 

* k = 8, 16: rows 4,8,12,... are also 0, ie. g 2 'd+4m = mod 2 l (/ < 4) 

* k = 32: no rows anymore! 

* k = 5: all rows are 0, idem at k = 25. 

* k = 20: row 4,8,12,16,20 are 0, because both are mod 4 and mod 5. 

These congruences are much less impressive than in the previous cases, due to the more 
complicated origin of the instantons. It is not even mathematically understood why 
these are integers and what exactly they count. Again we have not attempted to prove 
the congruences. 

5.2.1 Asymptotics 

In this case the asymptotic behaviour is much more tricky than in the previous two 
examples. The growth is indeed exponential, but logg^ has more than a simple loga- 
rithmic term and monomial terms. Indeed, subtracting the logarithmic term gives us a 
sequence on which the asymp^. trick works badly - as if other logarithmic terms were 
hiding. In fact, finding out the coefficient of the first log term is already quite tough, and 
differentiating (to get rid of log-terms) does not yield anything with only monomials. 
The second author's attempts to deal with the asymptotics of log qa through the asymp^ 
trick can be found in [TTj . 

6 Conclusion 

Though the congruences satisfied by the sequence v n of lines in P n are numerous and 
very interesting, the exponential asymptotic growth v n ~ n 2n rules out that the v n are 
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Fourier coefficients of any modular form on a subgroup of SX(2,Z) (whose coefficients 
typically grow like n 2k ~ l for weight 2k). It is quite gratifying to see that the asymptotic 
expansion in (|1U|) can be written out with as many exact terms as one wishes, since the 
coefficients are rational. 

Another sequence - that of numbers rid of degree d curves through 3c? — 1 points of 
the plane - has very similar behaviour, in terms of congruences as well as asymptotics 
(though the latter's coefficients will be irrational and not recognized). A last sequence 
- that of numbers of degree d rational curves on the quintic Calabi-Yau threefold - is 
much less enticing; its congruences are rather limited and its asymptotics are awkward: 
not simply one or two log-terms followed by mere monomial terms, but certainly log(log) 
terms or infinitely many log terms. 

It would be interesting to study more of typical sequences from enumerative geometry 
and see if there is an underlying pattern. Also, the question of how many of those 
sequences satisfy a recurrence relation is still open. Among the three sequences that we 
discussed, only that of nd (plane rational curves) obeys a known recurrence. 
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In this appendix we prove the alternative definition (j3J) and the asymptotic formulae (j3J) 
and (fTU|) for the numbers v n defined in (J2J). 



Proposition 1 Let G(x,y) be a homogeneous polynomial of degree 2n in two variables 
and P{x) a monic polynomial of degree n + 1 with distinct roots. Then the expression 



APPENDIX 



by Don Zagier 
Exact and asymptotic formulas for v. 



Exact formulas 




(12) 



P( a )=P(/3)=0 



is independent of P and equals the coefficient of x n y n in G(x,y). 



24 



Proof. By linearity it is enough to consider monomials G(x, y) = x r y s , r + s = In. 

(a)=0 p'( a ))Kl^P(f3)=0 P>(f3) ■ 



Then the expression (fT^j) factors as (Y^p( a )=o p^h))(^2p(/3)=o pkffd- by ^he residue 



theorem we have 

'x r dx\ „ (x T dx 

IXGS x=a I 

P(a)=0 ~ a6C 



= R (x[dx\ fx^dx^ 
2^ pr( a ) 2^* eS *=«{p( x )) * eSx =°°{p(x)J 



and this equals if < r < n and 1 if r = n since P is monic of degree n + 1. The 
proposition follows. □ 

Remark . The same proof shows that if G is homogeneous of degree m + n and P and Q 
are two monic polynomials of degrees m + 1 and n + 1 with distinct roots, then 

^ P'(a)Q'(P) 

P(a)=Q(/3)=0 V 7 ^ Kf J 

is independent of P and Q and is equal to the coefficient of x m y n in G(x,y). Yet more 
generally, and still with the same proof, if G is a homogeneous polynomial of degree 
ri\ + • • • + rife in k variables and Pi, . . . , P& monic polynomials of degree n% + 1, . . . , + 1 
with no multiple roots, then 

G(o;i, . . . , a k ) 



P'(ai) • • • PL(a k ) 

P 1 (ai)=...=P k (a k )=0 U 1J kK kJ 

is independent of all the Pj and is equal to the coefficient of x" 1 • • • x^ k in G(xi, . . . , 

In fact G need not even be homogeneous, but can be any polynomial in k variables of 

degree < n\ + V n k . 

Corollary 2 Let F(x, y) be a symmetric homogeneous polynomial of degree 2n — 2 in 
two variables and Wq, . . . , w n distinct complex numbers. Then the expression 

F(wi,Wj) 



£ — 

0<i<j<n PI {Wi-W k ){Wj -W k ) 



0<k<n 



is independent of u> , • • • , w n and equals the coefficient of x n in (1 — x)F(x, 1). 

Proof. This follows after a short calculation if we apply the proposition to G(x,y) = 
{x-y) 2 F(x,y),P{x) = Yl n i=0 {x-Wi). □ 

Corollary 2 immediately implies that the right hand side of equation (J2J) is indepen- 
dent of the (distinct) complex variables w ,...,w n and that (0) is equivalent to (jSJ)- 
The computational advantage is huge: formula (0) is very slow to compute, even for 
moderately large n, whereas © can be implemented in PARI in one line as 

v(n) = coeff (prod(j=0,2*n-3,2*n-3-j+j*x, 1-x) ,n-l) 

and takes < 2 seconds to compute v n up to n = 100 and 46 seconds up to n = 224. 
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We can rewrite (JBJ) in several other forms by using residue calculus. Setting D = 2n—3 
and making the substitution x = 1 — D/z, we find 



2n-3 



Res x=0 ((l-a;) J](2n-3-j+j 



j=0 



17 2 » R<*. : L) ( U t° {Z j <h). 



dx\ 
x) — 



z n+i ( z _ D y 

Since the residue of the integrand at infinity is zero, we can also write this as 

v n = - D 2n Res 2=0 



(13) 



(14) 



D 2n Res, =n ( -^*F^ ^- dz 



z n+i ( z _ D y 

while simply making the substitution z \— > D — z in (|14|) gives the similar expression 



(15) 



,z n (z- D) n+1 

and adding these two last expressions gives yet a third form: 

i„^^i„ / rii=o(^ _ j 



(16) 



v n = -D 2n+1 Res z=0 



z n+l U _ D )n+1 



dz ) . 



(17) 



Each of the formulae ([15 p -([17j ) expresses v n as the constant term at z = of the Laurent 
expansion of a rational function, e.g. (JT5j) says 



-l) n D 2n ■ coefficient of z n ~ l in ^^^j^" 1 "^ as : 0. 



;is) 



Substituting z = Du, we can write this as 



v n = {-l) n D 2 ■ coefficient of u n ~ l in ^-Du)(2-d^(d-i-Du) as m ^ 0; (19) 

from which we see again that D 2 \v n (Lemma QJ. By expanding (D — z) l ~ n by the 
binomial theorem, we also obtain closed formulae for v n ; for instance (fT%|) gives: 



n-l 



Vn = Y j {-i) n - 1 - m ( 2n ~ 2 ~ m \[y>^ 

m=0 



n — 1 



D 

m 



(20) 



where 



D 



first kind. 



the coefficient of z m in z(z + 1) • • • (z + D — 1), is a Stirling number of the 



Asymptotics 

To obtain the asymptotic expansion of v n , we write the residue in (II 3|) as ^ J| a .| =1 an d 
we make the substitution x = (1 + it)/(l — it) to obtain, after a short calculation, 



-j-j- ^D 2 + r 2 t 2 \ t 2 dt 



o /*oo 



t 2 dt 

i + t 2 ) 2 ' 



(21) 
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where </>d(£) denotes the rational function 

1 + r 2 D~H 2 



Mt) = JI 



1 + t 2 

r=l,3,...,D 



It is easy to see that 0r>(O) = 1 and (pD(t) < e~ cDt2 for some absolute constant c > 
(a much more precise formula will be given in a moment), so the main contribution to 
the integral comes from small t. For t small and D large we have (uniformly in both 
variables) 



3=1 J r=l,3,...,C 



f 2 ' 



D 1 \ 2 (D 1 2 \ , ( D 1 4 / 1 \\ fi 

-(f-^^^Q)^(-n^^Q)' 10 

■ M\ .1, . ( D . ^1^ . 14 , (D , „(\ 



V13 3D 
and hence 



°(^)) iI2+ (-IS + (^)) fl4+ (T7 + (^)) fle + 0(fl(1 " 



2 

x , , a; 



■ x 2 /D) 2 y/D' 

-x 2 /3 



e 



x 2 + (— - 2x 4 )/T 1 + f- ^— + 3x 6 + — 

L V 5 / V 50 35 3 / 

/x 14 12x 12 314x 10 59x 8 fi \ , 

+ 1 x 6 )D~ 3 H 

V750 175 315 15 / 

a; 30 Hx 28 3 55x 10 2x 8 \ 7 

+ 1 1 1 lD~ 7 + 0(D~ 

V 393750000 19687500 162 45 / v 



Substituting this expansion (with the 34 omitted terms included) into equation (|2ip 
with t replaced by xj\fD and using the standard evaluation 



n! V4/ 



we obtain 



„„ = M D °-V*(l - + gggfl-^ - + 25225773 4 

V 7T V 4 160 3200 358400 

10092025737^ = 2271842858513^ fi 4442983688169 ^ , „, 

35840000 2007040000 1146880000 v 

This asymptotic formula can of course be written in many other ways, e.g.: 

87261 



or 



9 


111 


9999 


8n 


640n 2 


25600n 3 


15 


1689 


79281 


8^ 


640n 2 


f 25600n 3 
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5734400™ 4 
19691853 



e-\ -(2n) 2 "-^ 2 1 + — + — — + + — — — r + 



v n , n 11 141 9973 59673 

log - — - = 2n - 4 \ogn + C H 1 H H r H 

6 (2n)! 6 6n 160n 2 28800n 3 179200n 4 

with C = — 3 — log7r — | log |. Of course, more terms could be obtained in any of these 
expansions if desired. 
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